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Abstract 

We report the discovery of a short-duration microlensing candidate in the Northern field of the POINT-AGAPE 
pixel lensing survey towards M31 (a = 00'M2'"51.42'' and 6 = +41°23'53.7"). The Full-Width at Half-Maximum 
(FWHM) timescale is very short (tf„hm — 1-8 days) and the source star has almost certainly been identified 
on HST archival images, allowing to infer an Einstein crossing time : Ie ~ 10.4 days and a magnification at 
maximum : Ams,^ ~ 18. 



1 Introduction 



In 1986, B. Paczyhski suggested to search for 
MACHOs using gravitational microlensing on back- 
ground stars. The possibility of detecting such events in 
M31 was independently suggested by A. Crotts and 
P. Baillon et al. [Q. The advantages of targeting such a 
large external galaxy are : 

• the number of stars that can act as possible sources is 
enormous. 

• this gives access to a different direction than 
Magellanic clouds to study the halo of the Milky Way. 

• the high inclinaison of M31's disk causes an asym- 2 Pixel lensing 

metry in the observed rate of microlensing by lenses in 
a spheroidal halo, as shown in figure This gives an 
unambigous signature of the halo lenses 0] . 



called "Pixel method" (see section 2.1) and identi- 
fied several candidate microlensing events. In par- 



ticular two short-duration candidates was reported : 
AGAPE-Zl § and PA-Nl g. The Columbia- VATT 
group 1^ developped independently a method based on 
differential image analysis. As in most cases we do not 
have acces to the magnitude of the sources, there is a 
parameter degeneracy in the theoretical "Paczyhski fit" 
and the way to compute the optical depth is more com- 



plicated than for resolved stars (see section 2.2) 



2.1 Pixel method 
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Figure 1: The high inclinaison of M31's disk causes an 
asymmetry in the observed rate of microlensing events 
by lenses in a spheroidal halo H. 



The main difficulty of such observations is that 
the sources are not resolved (nevertheless they might 
be resolved when they are lensed). The AGAPE 
collaboration 0| has developped an analysis method 



The Pixel method was suggested by P. Baillon et al. ^ 
and A. Gould §. The light-curve of every pixel is stud- 
ied to detect significant variations. Due to the seeing, 
the flux of a star is spread over several pixels, as shown 
on figure ||. Therefore several neighbouring pixels are 
simultaneously affected by a single variation. To get 
the best signal to noise ratio it is necessary to find the 
central pixel of an event. The event is detected if the 
increase of the flux in the central pixel is typically bigger 
than the noise : 



^^central 
pixel 



[A - 1] X / X > 



where A is the magnification, / the seeing fraction and 
the fiux of the source. As only large amplifications 
are detected, the detection efficiency is weak but this is 
compensated by the huge number of monitoring sources. 



To appear in the proceedings of the Dark Univers conference, 

2-5 april 2001, Baltimore (http://ntweb.stsci.edu/SD/DarkUniverse). ed. Mario Livio. 
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Figure 2: The pixel method consists in the survey of the 
flux in each pixel. An event is detected is the increase 
of the flux in the central pixel is typically bigger than 

the noise : A^^entral ^ cr- 
pixel 
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Figure 3: In pixel lensing one can usually determine 
= X [A(tmax) - 1] - *bkg but not 



*bkg and $ 

ma; 

and so not Ie- 



2.2 Optical depth in pixel lensing 

The flux increase in a microlensing event is given by (for 
pointlike lens and pointlike source) : 



A$(<) 
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$^ is the source flux in absence of lensing. 
u{t) = 6{t)/6E; 0{t) is the angular separation of lens and 
source, and 6e is the angular Einstein radius : 



9e = 



AGM s - I 



si 



where M is the lens mass and s and I the distances of 
the source and the lens from the observer. The rela- 
tive motion of the lens and the source is usually well 
approximated by a uniform motion : 



P tE 



tmax is the time of the maximum magniflcation, /3 the 
impact parameter in units of the Einstein radius and 
tE is the Einstein time. In the limit of large magnifica- 
tions at maximum, the microlensing curve assumes the 
degenerated form : 
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The three parameters f3 and tE have reduced to 
two: "$H./;9" and "/Jts". A fit to the data can deter- 
mine neither the star flux nor the Einstein time tE 
separately from the unknown impact parameter jS. This 
degeneracy is illustrated in figure ^. 

Usually, the estimated optical depth is defined as : 



2NsT ^ e, 

detected 
events 



where T is the estimated optical depth, the number 
of observed sources, T the duration of observations, i^;. 
the Einstein time of event i and the detection cfli- 
ciency for event i. As it is shown above, it is usually 
not possible to determine tE in pixel lensing. Remark- 
ably, however, it is it still possible to evaluate the optical 
depth. P.Gondolo 1^ shown that the estimated optical 
depth can be written : 



0.72 
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where 3*1^'^ is the detection threshold and Fgg a func- 
tion of ^Yirn^ which depends slightly on the luminosity 
function of sources. 



3 POINT-AGAPE collaboration 

The POINT-AGAPE collaboration uses the Wide Field 
Camera (WFC) on the 2.5 m Isaac Newton Telescope 
(INT) at La Palma, Canarie Islands . It carries out 
a pixel-lensing survey, monitoring two flelds of 0.3 deg^ 
each, located North and South of the M31 center. The 
pixel size is 0.33". The survey has the potential to map 
the global distribution of the microlensing events in M31 
and to determine any large scale gradient. Three years 
of data are planned (from August 1999 to January 2002). 

3.1 Data analysis 

For the moment we restrict the analysis to the Northern 
field in 1999, as shown figure |[ The observations are 
spread over 36 epochs between August and December. 
The exposures are in two bands : 36 nights in Sloan r' 
and 26 in Sloan g'. The exposure time is typically be- 
tween 5 and 10 minutes per night and per band. 

After bias subtraction and fiat-fielding, each image is 
geometrically and photometrically aligned relative to a 
reference image (14 August 1999) which was chosen be- 
cause it has a long exposure time, typical seeing (1.6") 
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Figure 4: The two fields observed by POINT-AGAPE 
collaboration 



and i? = 20.87 ±0.13 V - R = 1.2 ±0.22 which is 
very red. On figure (panel a) one can see that, after 
t ~ 120 days (i.e in December 1999), a second bump ap- 
pears on the superpixel light curve. Using image sub- 
traction, one can show (panels d and e) that this second 
bump is not due to the same source as the first one. Dis- 
tance between the two sources is about 1.1". December 
points were excluded from the fit. 

The PA-Nl position lies within a serie of five Hubble 
Space Telescope (HST) WFPC2 archival images taken 
in July 1996, three with F814W filter (~ I) and two 
with F606W filter (~ V). On these images we found a 
resolved source close to the la error box and no other 
within the 3cr error box. The prior probability to find 
such a star so close to the predicted position is only 3%. 
The color of this resolved star (F = 24.51 ± 0.12 and 
1^ — / ~ 2.10 ± 0.16) is perfectly compatible with the 
color of PA-Nl. Moreover figure @ shows that the source 



and little contamination from the Moon. The light 
curves are made for (2.3" x 2.3")-superpixcls to cope 
with seeing variationq^. 

We use a simple set of selection criteria designed to 
isolate high signal-to-noise events. Detection of events 
is made in the Sloan r' band which has better sampling 
and lower sky background variability. The Sloan g' band 
is used to check achromaticity. 

3.2 Preliminary results 

In 1999 data and in Northern field, we have already 
found about 300 microlensing candidates (compatible 
with a Paczynski fit) which will have to be validated 
using 2000 and 2001 data. Magnitudes at maximum of 
amplification are in the range : 21 < Rmax ^ 24.5, where 
R is the standard Johnson-Cousins band. 

One of our candidates, called PA-Nl, is par- 
ticularly interesting with ifwhm < 2 days and 
A$max ~ 18 ADU.s~\ corresponding to i? 20.8. 
This candidate, presented in [^, is the subject of the 
next section. 

4 PA-Nl microlensing candidate 

Figure || (panels a and b) shows the lightcurves in 
Sloan r' and Sloan g' of this candidate together with 
the theoretical fit derived below. Using the Aladin Sky 
Atlas O, we found that PA-Nl has J2000 position : 
a = 00M2'"51.42" and S = +41°23'53.7". That is, it lies 
at 7'52" from the center of M31. On images close to 
the maximum amplification, PA-Nl is bright enough to 
perform its photometry. We found : V = 22.00 ±0.17 

^This size is set by the worst seeing (~2.1"). 
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Figure 6: Superposition of the Hipparcos (points) and 
resolved stars in HST images (shaded area) color- 
magnitude diagrams. The Hipparcos stars was moved 
on the position on M31. The cross hairs show the 
position of the HST source {V = 24.51 ± 0.12 and 
T/- 7 = 2. 10 ±0.16). 



star is either the star visible on HST images or must be 
on the main sequence and so magnified by Amax ?J 10^- 
The latter possibility is extremely unlikely. Hence a firm 
prediction of the microlensing interpretation is that the 
source is almost certainly identified to be the red giant 
visible on HST images. Using this information, we can 
eliminate the degeneracy between fit parameters. The 
parameters are shown in tabic pi. 
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Figure 5: Panels a and b : flux in r' and g' against time in days; days = J-2451392.5 where J is the Juhan 
date, i.e: data are spread over August and December 1999. Panel c : zoom centered on the event showing the 
variation of the ratio of the flux change in the two passband with time. The vertical hues are centred on tmax and 
are separated by 0.9 days, i.e. 1/2 FWHM. Panels d and e : positions (using image differencing) of PA-Nl and 
the December bump. One can see that the two sources arc clearly separated by > 3 pixels, i.e ~ 1.1" (pixel size 
is 0.33"). However the December source is close enough to PA-Nl to be visible on the supcrpixel lightcurve. 



^max 

(days) 


/3 


tE 

(days) 


(ADU/s) 


(ADU/s) 


'i'bkg.r' 

(ADU/s) 


*&bkg,g' 

(ADU/s) 


13.87 


0.056 


10.4 


1.02 


0.28 


397.65 


209.4 



Table 1: Fit parameters assuming the amplified star is the red giant visible on HST images. 
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